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We study self-consistently the finite T broadening of tfie p-meson and its implications for dilepton 
emission in heavy ion-collisions. For this purpose finite width effects at finite temperature due to 
the p-TT interaction are investigated in a self-consistent ^-functional approach. The temperature 
dependence of the p-meson and pion spectral functions and self-energies is discussed. The spectral 
functions show considerable broadening in comparison with a perturbative calculation on the one- 
loop level. Using these spectral functions, we make a comparison to dilepton emission data from 
the CERES NA49 collaboration employing a parametrized fireball evolution model of collision. We 
demonstrate that these non-perturbative finite width effects are in-medium modifications relevant 
to the understanding of the enhancement of the low invariant mass spectrum of dileptons emitted 
in A-A collisions. 



I. INTRODUCTION 



In recent years many theoretical and experimental research has been undertaken to investigate the physics of ultra- 
relativistic heavy-ion collisions (URHIC). The important goal here is to gain a better understanding of the properties 
of extremely hot and dense nuclear matter. Through variation of the experimental parameters - e.g. collision energy 
and impact parameter - different regimes of temperature and net-baryon density of the nuclear system are explored. 
The ultimate goal is to understand the phase structure of quantum chromodynamics (QCD), the theory of strong 
interactions. 

In order to extract reliable information about the properties of nuclear matter under these extreme conditions (and 
especially to decide if a phase transition from the phase of (confined) hadronic matter to the (partonic) quark-gluon- 
plasma (QGP) has taken place) one has to do a careful and comprehensive analysis of the experimentally accessible 
observables. Dileptons (e^e~ and fi^ fi~ pairs) and real photons are interesting probes in this context. Since they 
are emitted during all collision phases and leave the interaction region without substantial rescattering due to the 
smallness of the electromagnetic coupling. Hence, they reflect the vifhole space-time evolution of the collision process, 
in contrast to hadronic probes which predominantly reflect the conditions at kinetic decoupling. 

To some degree, measurements made so far complement each other: while photons have been measured for compar- 
ative large transverse momentum scales of 1.5 GeV and above, the invariant mass spectrum of dielectrons is rather 
well known below 1 GeV. Therefore, the measured photons naturally carry information about the early, hard partonic 
stages of the collision process whereas the scale set by the dilepton measurement dominantly reflect the later, soft 
hadronic stages. 

While there is little reason to assume that low mass dilepton physics would be dominated by emission from a 
QGP, understanding the late hadronic phase of the evolution is nevertheless of crucial importance for forming a 
comprehensive picture of the whole collision process. The most intriguing physics question here is the possibility of 
in-medium modifications of hadron properties. The CERES and IIELIOS-3 collaborations demonstrated that central 
nucleus-nucleus (A-A) collisions exhibit a strong enhancement of low-mass dileptons in the range between 200 MeV 
and 600 MeV as cornpared to the yield one expects from the corresponding rates of proton-nucleus (p-A) and proton- 
proton collisions Q, ^ IS which may hint at such modifications. For an overview of the various mechanisms 
brought forward to explain this enhancement see e.g. 01 • 

It is often argued that a large broadening of the p-meson spectral function cannot be expected by a purely mesonic 
medium. This assumption is based on perturbative calculations of finite width effects of pions on the p-meson at finite 
temperature However, these investigations consider the width of the p-meson caused by stable particles only. In 
most of the theoretical investigations the damping width due to collisional broadening is ignored or investigated in 
an extended perturbative frameworkj^. Although investigation of collisional broadening effects due to the interaction 
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of vector mesons with thermal mesons in a kinetic theory framework have been carried out, they consider possible 
scattering processes of on-shell mesons only [lOl ITU . In contrast we extend these studies in order to include also 
off-shell contributions. 

In the present paper we study these broadening effects and demonstrate that even the n - p interaction alone is 
able to induce substantial broadening of the p-mesons's spectral function due to the in-medium damping width of the 
pion. 

For the purpose of studying this mechanism, we employ a many body rcsummation which is known as the 
functional approach. In this approach we derive Dyson Schwinger equations for the propagators of the p-meson and 
the pion. Our calculation arc self-consistent in the sense that the in-medium broadening of the p-meson and the pions 
is taken taken into account by this scheme in its reciprocal dependence. 

We concentrate on the phenomenology of coUisional broadening and leave the effect of mass shifts relative to the 
vaccum mass of the particles (especially problems of renormalization) aside by taking into account only the imaginary 
parts of the self-energies. 

Since our primary aim in this study is to investigate how the tt - p interaction, giving raise to the dominant decay 
channel of the p-meson in the vacuum, is modified at finite temperature, we disregard other effects in our model (like 
mass shifts, scattering off by baryons and effects of chiral symmetry restoration) focus sing on coUisional broadening 
due to that interaction only. For approaches considering these other effects, see e.g. [7l ll2l[T3j . 

We compare the results of this self-consistent calculation with a perturbative analysis and show that the self- 
consistent treatment of this interaction leads to considerable coUisional in-medium broadening at finite temperature. 

In order to demonstrate that the broadening effects are quantitatively important for the understanding of the 
observed low mass dilepton enhancement as seen in experiments, we compare the model predicitons with the exper- 
imental results obtained by the CERES-coUaboration 0, 0, Q . For this purpose, we fold the predicted rate with a 
parametrized fireball evolution which has been shown to successfully reproduce a number of different observables, 
among them charmonium suppression, direct photon emission, hadronic transverse momentum spectra and Hanbury 
Brown Twiss correlation measurements [T^ . 

The paper is organized as follows: We first introduce the framework for the description of the n-p interaction. Then, 
the self-consistent resummation scheme is described and the approximations of the Schwinger-Dyson equations for 
the p-meson and the pion are derived where we employ the Stiickelberg method to describe the p-mesons as massive 
vector mesons. Special care has to be taken to propagate only four-dimensionally transverse and therefore physical 
modes of the vector-mesons in order not to violate current conservation. Two different solutions of these difficulties 
are discussed. In section Hvl we present the solution of the self-consistent equations for the spectral functions of the 
pions and the p"-meson and discuss the quantitative and qualitative features of the in-medium modifications which 
are important in the hadronic phase. Section IV Bl contains a brief review of the model description of the QGP phase, 
section IV Dl deals with non-thermal contributions to the dilepton yield, section IV Al provides a short introduction 
to the parametrization of the fireball evolution. Finally, in section I VII we compare the model calculation with the 
dilepton spectra as obtained by the CERES collaboration and discuss these results. We end with a conclusion and an 
outlook. 



II. THE p-MESON tt INTERACTION MODEL 

We restrict our self-consistent study of the in-medium properties of the p-meson to the p-ir interaction, emphasizing 
self-consistent finite pion width effects on the p-nieson. The framework of these considerations is a model inspired by 
vector meson dominance (VMD) 0,^^ for the interaction of neutral p-mesons with the charged pions. 

In order to describe the propagation of the massive vector mesons in the appropiate quantum field theoretical 
framework, in order to propagate their physical modes, we treat the p-meson in the Stiickelberg formalism [13. (For 
a detailed review of this formalism see also The suitability of this formalism for calculations in ^-functional 

approaches was first discussed in |l9l |. In the Stiickelberg framework one introduces an additional scalar field, the 
so called Stiickelberg ghost (j>. The free Langrangian density of the p-meson field and the Stiickelberg ghost cj) is 
given by: 

1 1 

-Cstucck = -^^M-^'"' + — ^M^"" + ^^d^^'^^d^'^) + mc^d^A'^ , (I) 

where = dfj,Af^ — d,yAf^ is the corresponding field strength tensor. Pauli showed first [23| that this treatment 
of massive vector mesons leads to a symmetry of the action (up to a total divergence) under the following (local) 
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transformations of the fields:. 

Af, ^ A'^ = Af, + SAf, = Af, + df,X , 
(f) ~* (f)' = (f) + 5(t) = 4) + m\ . (2) 

This gauge symmetry has to be taken into account when quantizing the theory. Applying the standard Fadeev 
Popov quantization scheme and introducing a gauge fixing term in the i?^-gauge 

g{A,<i>)^^^A^^ + ^,m^ (3) 

one obtains after quantization the following effective free Lagrangian density including Stiickelberg ghosts (p and 
Fadeev Popov ghosts -q: 

= + (a^ - -d^'A - 4 + ^rncl^f 

4 2 \ my 2^ 

-77* + ^m^) 77 . (4) 

Expanding the contributions in the effective free Lagrangian density and partial integration shows that in the 
i?j-gauge the fields A^, cj) and 77 all decouple from each other. 

There are many advantages of this description of massive vector mesons. The theory is symmetric under the appro- 
priate Becchi-Rouet-Stora-Tyutin (BRST) - transformations |^ l^^, and leads to a positive definite Hamiltonian 
for the physical propagating mo des jlSl. Moreover this leads to proof of the unitarity and can be used to demonstrate 
(perturbative) renormalizability |l8ll24j . 

The neutral p-mesons are in a VMD-inspired model decisive for the dilepton production, which is given in the hot 
p — TT medium by the decay process p*^ — > 7* ^ e"^ + e~ via an intermediate virtual photon 7*. We now include the 
interaction of the neutral p-mcsons with the charged pions. 

The pion is coupled minimally to the vector meson field and to the electromagnetic photon field Af^_^. Therefore 
the part of the Lagrangian density including these interactions is given by: 



— miTT*TT . 



(5) 



Additionally, we introduce an interaction term coupling the photon to the p-meson vector field; 



where !F^'^ and i^'"' are the field strength tensors of the rho-meson and photon, respectively. 
The total Lagrangian density ^ is now the sum of the three contributions discussed above: 

The electromagnet current created by the mesons is given by: 

4 = %d^Fi^ + -m'A, (8) 
9 9 

and is in order e/ g proportional to the p field. This is the manifestation of VMD. 

The coupling of the photons-field to dileptons can be calculated in quantum electrodynamics (QED). This is essential 
for describing the production of dileptons from the hot hadronic medium in VMD, see appendix IbI 



^ It should be noted, that if one neglects all terms including the field, this Lagrangian for the pions and the p^-meson is formally 
identical to that one for scalar quantum electrodynamics with a massive photon. The pions as pseudoscalar have taken the role of the 
scalar particle and the p^-meson has taken the formal role of a massive photon. 
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III. SELF-CONSISTENT RESUMMATION SCHEME 

We employ a self-consistent many body resummation scheme which has been first introduced by Baym and is 
known as the <f>-functional formalism to model self-consistently finite width effects on the p-meson by the pion due to 
collisional broadening, ft has been extended to the case of relativistic quantum field theory by Cornwall, Jackiw, and 
Tomboulis 

This functional method is based on the resummation of the partition sum. In this scheme one obtains Dyson 
Schwinger equations for the full propagators of the p-meson and the pions with the following form: 

S;.' (K) = {!<) + n^, (K) , (9a) 
r-\K) = p-^iK)+U{K) , (9b) 

where S and V represent the full propagators of the p-mcson and the pions, respectively. S and P are the free 
propagators of the p- meson and the pions: 

S;,\K) = [K^ - ir?)g^, + ^-^K^^K, , (10a) 

P'^{K) = -ml . (10b) 

The self-energies 11^'' and 11 are the self-energies of the p-meson and the pions, respectively, which are given 
as functional derivatives of the sum of the two-particle irreducible diagrams V2 of the theory in this resummation 
formalism: 



n(P) . 2^^^[^'^] 



n^.(P) = 2 



SP{P) 
6V2[S, P] 



SSt^'iP) 



and (lla-) 

S=S,P=V 

(lib) 

S=S.P=V 



Taking all two-particle irreducible diagrams into account corresponds to solving the full quantum field theory. 
A systematic approximation scheme is obtained by taking into account certain classes of two-particle irreducible 
diagrams. In our analysis we consider the following approximation of V2: 

V2 = -g' ^ l^{2P + Lr{2P + LYS^,{L)P{P)P{L + P) , (12) 

which corresponds to the diagram on the left of fig. ^ Our notation is J^, /(fc) = T^^^_^ J (i'^k/(27r)'^ /(27rmr, k). 
This scheme does not include any resummation corrections of the vertices, so in V2 one takes the vertices of the tt — p 
interaction stemming from the Lagrangian density JT)) on the tree-level. 

Calculating the self-energies for example for the p-mesons according to eqn. Hll|l corresponds on the diagrammatic 
level to cutting the p-propagator in the diagram of V2. In an analogous way one obtains the self-energy of the pion 
(see fig. ^ . The scheme is self-consistent. Because every propagator appearing in the self-energy contributions is a 
full-propagator, i.e. in a pcrturbativc language the propagators in the loop diagrams are fully resummed on that loop 
level. 

The self-energies of the pions and the p°-meson are therefore given by (|11|) : 

n(P) - 2^1^ - -4g2 J (2P - Qr{2P - QrS^,{Q)V{P - Q) 

(13) 

and 

np.(P) - ^ J^\p^ = -2.9' ^(2Q - P)m(20 - P)uV{Q)V{P - Q) . (14) 

It is obvious that these self-energies themselves functionally depend on the propagators which enter the Dyson- 
Schwinger equations. 




FIG. 1: On the left the diagram representing the approximation for V2 is depicted. Via functional variation one obtains the 
self-energy contributions to the pion and p-meson self-energies as on the ride hand side. The wiggled line represent the p-meson 
propagator and the solid line the pion propagator. The diagrams have to be resummed self-consistently. 



As we are interested in the dilepton emission from a medium in thermal equihbrium, we need to calculate the 
retarded propagator of the p-meson in thermal field theory, see e.g. That quantity directly enters the calculation 
of the dilepton production rates. For that purpose we recall that the thermal average (A) of an operator A is defined 

by 



1 



Zif3) 



Tr(ic 



(15) 



where H is the Hamiltonian, (3 = 1/T is the inverse temperature and Z(^f3) is the partition sum. 
propagators of the p-meson and the pions are defined as follows: 



The retarded 



= -^{0{t)[A'^{t),A^O)])0 , 
They can be expressed in terms of the spectral functions of the p-meson and the pions: 



(16) 



5r(fco,k) 

^R(fco,k) 



— 00 
00 



2tt fco — Po + ity 
dp'o p(Po,l<) 
2tt ko - Pq + ir/ 



(17) 



where an analytical continuation from Matsubara frequencies to real energies has been performed. 

The vector meson propagator is a Lorentz-tensor. Lorentz invariance is broken because the 4- vector of the medium 
defines a preferred frame of reference and the tensorial densities (as spectral functions, self-energies, propagators) 
become functions of energy and of momentum. To identify the physically propagating modes of the vector meson, 
we project its propagators, self-energies and spectral function onto the four-dimensionally transverse and longitudinal 
subspaces. A detailed overview of the conventions of the projection formalism is given in appendix [XI . 

Every tensorial quantity can be decomposed into the four-dimensionally transverse components (proportional to 
the projectors A and B) and longitudinal components (proportional to the projectors C, E), see eqn. (|A9|) . 
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In general the spectral function components can be expressed via: 

1 /•= 
PA(fco,k) = -2Ini-A,,,(fco,k)7^«'"'(fco,k) = -Im / 
^ J —I 



-A nf^" 



PB(fco,k) 
Pc(fco,k) 



-2ImB^^(fco,k)i:>^^''(fco,k) = -2Im 

B^,(fco,k)p'^''(fco,k) , 
1 



' dk'o A^,(fco,k)p^"-(fc(),k) 
2 27r fco — fco + 

dfc^ B^,(fco,k)p^-(fc^,k) 
27r fco — /cq + i?7 



-2Ini 
1 



9Z,2jL2 ^^ 



C^,(fco,k)i?«^''(A:o,k) 



Im 



dfc^, C^.(fco,k)p^''(fc^,k) 



27r 



1 



91.2 jL2 



C^,(fco,k)p^''(fco,k) , 



/9£;(fco,k) = 



dfc^, E^,(fco,k)p^''(fc^,k) 



27r 



f 

-2IniE^,(fco,k)Z?«^'^(fco,k) = -2Ini / 

J — ; 

= E^,(fco,k)p^''(fco,k) . 
The retarded propagator can be rewritten in terms of these spectral fmiction's components: 

°' 7-00 27r ka - k'„ + iri 

Inverting the resummed propagator (see eqn. I9a|) of the p-meson and employing eqns. (|18|l leads to: 



PA 

PB 

PC 



-2Im-A^,(po,p)5'"'(po,p) 



-2Im 



p2 



2ImB^,(po,p)5^"(po,p) = -2Im [b-\P)] , 
= -2Im— i^C^,(po,p)5^''bo,p) 



-2^Im 



-m 



P2 _ ^ ^jj^ 



PE = -2ImE^,(po,p)5^"bo,p) 



-2CIm 



b{P)- 



where the function is given by 



h{P) 



p2 



Bh 



pgp2n2 



P2 _ ^^2 ^ 



(18) 



(19) 



(20) 



(21) 



and the projected components B^, B;,, Be, and Be (see appendix of the self-energy B^^^ are introduced. The 
quantity ^ is the gauge parameter appearing in the Stiickelberg method in the i?^-gauge (see eqn. (|2J). 

In spite of the truncation of the resummation (e.g. limiting the sum of the two-particle irreducible diagrams V2 to a 
subset of all possible digagrams) conservation of conserved currents stemming from global symmetries as long as they 
are realised as linear representations of the conserved group detailed balance and unitarity as well as thcrmodynamical 
consistence. This is guaranteed in ^-functional approximation schemes, as shown by Baym Bowever, in these 
schemes the vertex corrections necessary to guarantee Ward Takahashi identies for the propagator are not taken into 
account consistently. This shortcoming can lead to serious problems of the self-consistent treatment of vector particles 
on the propagator level [T9| . 

This problem has to be addressed in our scheme as well. In j23| it has been shown that in order to be consistent 
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with the corresponding Ward Takahashi identities the following relation: ^: 

P'S'l^'"'' = ^T^T^2 (22) 

for the (fully resumnied) propagator of the massive vector meson S^n has to be fulfilled in the Stiickelberg theory. From 
that condition one can conclude that in Stiickelberg theories (as in QED) the self-energy tensor of the vector particle 
has to be four-dimensionally transverse which corresponds to vanishing He and lie components of the self-energy 
tensor. 

This relation on the correlator level is spoiled in the employed ^-functional approximation, since only a partial 
resummation is performed. This fact is well known, not only for Stiickelberg theories |28l |. but also from Dyson- 
Schwinger approaches to the description of QCD phenomena [23 and stems from a violation of the gauge symmetry. 

One may think that these fundamental considerations of the violation of Ward Takahashi identies don't matter from 
a phenomenological point of view, because we are interested only in an apropriate description of the p-meson studying 
it's importance for the understanding of dilepton spectra and not about intrinsically complicated considerations on 
the violation of local symmetries in gauge theories. But we calculate the production of dileptons in a VMD approach 
where the p-meson couples directly to the virtual photon, these questions have to be addressed in detail. In part IbI 
of the appendix the equation for the rate of dilepton production is given. Eqn. HIj4|l shows that the current-current 
correlator is directly proportional to the retarded propagator of the /o-meson. Current conservation in QED which is 
essential for an appropriate description of the dilepton production tells that the spectral-function of the p-meson has 
to be four-dimensionally transverse: P^j.p^'^ = 0. 

To guarantee conservation of the hadronic electromagnetic current, we work in an appropriate gauge: the -R^-gauge 
taking ^ ^. This leads to the following modification of the eqns.: 

PA - -2ImiA^,(po,p)5^''(po,p) = -2Im (P2_^2^n,)"' , 

PB = -2ImB^,(po,p)5^"(po,p) = -2Im (P2_^2^nb)"' , 

PC = -2Im C^^{po,p)S'""{pa,p) = , 

PE = -2ImE^,(po,p)5'^''(po,P) = . (23) 

The components pc and do vanish and the transversality of the spectral-function P^Pf^, = and therefore current 
conservation is guaranteed. In that way we solve the most serious problem, mainly former violation of electromagnetic 
current conservation. On the level of the p-meson self-energy we still have non-vanishing 11^ and Ilg components, but 
this does not lead to problems calculating the dilepton rates. In the following we refer to this technique to guarantee 
the transversality of the spectral functions as " approach F . 

In |28l | the authors use a different ansatz by forcing the self-energy of the p-meson to be explicitly transverse via 
the eqns. lfn2)l (with X = n), which are violated in ^-functional approximation schemes. In this ansatz one argues 
that the temporal components of the self-energy tensor are tied to the conseration of charge and therefore involve an 
infinite relaxation time of the conserved quantity in the full theory. The self-consistent results are different because 
they reflect the damping time of the propagators in the loop. This behaviour has been studied in detail in |3ll |. where 
it has been shown that current conservation could be restored via a Bethe-Salpeter ladder resummation. The spatial 
components of the self-energy tensor are in general be less affected by the resummation in the case that the relaxation 
time for the transverse current-current correlator is comparable to the damping time of the propagators appearing in 
the loop, whereas the time components are expected to suffer significant changes |28| . This physical reasoning leads to 
enforcing the transversality of the self-energy tensor via the first two eqns. of ljA12p (with X = 11), which give the p- 
mcsons four-dimensionally transverse self-energy components (Ha , lib) as functions of 11* (P) = i [5^^ — fp^ (P) 
and n^(P) ^ p,W^P)pj: 

n^(P) = -n*(P) , (24a) 

n'^(P) = -^Il^(P) . (24b) 
Po 



2 This relation is somewhat similar to those for QED, e.g.P'^5^™ = ^P„/P'^ and QCD, e.g. Pi^P^S^^T^ = f- 

^ That is similar to taking the landau gauge in Dyson-Schwinger approaches to QCD at finite temperature, see e.g. |30 
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These equations are equivalent to eqn. (4) in |[28j up to a different sign convention. In the following we refer to this 
technique to guarantee the transversality of the self-energy tensor first employed in |28l| as " approach IF . 

For the sake of simplicity we are concentrating on damping width effects of the particles and avoid renormalization 
which would require a temperature independent self-consistent subtraction of infinities |32l | and don't take into account 
any changes in the real part of the self-energies. This ansatz to disregards any mass shift of the particles is analogous 
to the one in [28l |. 

Given the fact that we are interested in the phenomenology of the broadening due to finite width of the pions this 
approximation of neglecting mass shifts (i. e. disregarding renormalization), can be justified. In the n-p interaction 
has been studied perturbatively. The value of the mass shifts (relative to vacuum) predicted from this perturbative 
study are comparable with the mass shifts one can expect in any many body scheme. (This is different regarding the 
coUisional broadening, since in the perturbative calculation the width of the pion giving raise to the p meson's width 
is zero.) Especially from Fig. 6, 7, 8 and. Fig. 9 in Q one expects physical mass shifts from the 7r-p-interaction only 
in the range of about 10% of the vacuum mass of the p-meson even for comparably high temperatures. 

In that way, the equations for the non-vanishing components of the spectral function are: 



PA 

PB 

PC 

PE 



-2ImnQ 



(p2 _„j2)2 + (Ininb)2 






(25) 



To close this self-consistent set of equations, we have to express the self-energies as functionals of the spectral 
function of the pion and the spectral function's components pa and pb of the p-meson. We use the decomposition of 
the spectral function of the p-meson: 



p^^iP) = A^,Pa{P)+B^,pb{P) 



Sij + PiPj I Pa 



p2 



PtP] ]ph{P) , 



where pi is the normalized i.th component of the momentum vector pi = Pi/p and the following contractions 



MP-QTi^P-QT = 2iP^ + Ll)-Ll 

/ u2 r2\2 



^^{2P - Qr{2P - Qr 



(p2 



5,,{2p^qy{2p^qy = 2{p' + ll)~l 



2 

1 ; 



^-l^{2p-qr{2p-qY 



(26) 



where L\ ^ Q and L2 = P — Q. The quantities /^,^2j ^^'^ ^'"'^ ^^le squared three-momenta, /lo and Z20 are the 
energy components of Li and L2, respectively. We obtain the following equation of part of the self-energy of the pion: 



01 2 I 72\ ,2 {P" %Y 
2(P +'2) - '1 72 



Pa 



2^2 



- [2(P2 + Ll) -Ll + 2(/ + ll) -ll- '^II^'^ 

( 2 _ j2\2 > 
~^^^T^]Pb{ll0jl)jp{l20j2) 

X e(|/i-/2| <P<h+l2) 

X [1 + film) + fiho)] S{po - Im - ho) ■ 



(27) 



To determine the imaginary four-dimensionally transverse components of the self-energy Imlla and Imllf,, we employ 
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eqns. (jAll|) fwith X = 11) and use the following contactions: 



{2q-pf{2q-pf = iho-hof , 



{2q-pf{2q-p\^ = - 



{ho — ho){li — I2) 



We need to calculate Im ((5*^11^), Im (p^p^II*-'), Im (ll™) and Im (n*°pi). These are given by: 



(28) 



and 



Im {5'm,_ 



Im {p^p,\r') 



Im fn' 



00 \ 



p{2n) 



2ti 



d^20 

2tt 



X [2(2^ + ll) - /] p(/io, ;i)p(/2o, h)5 (po - ho - ho) 
X eilh-hl <p<li + l2) 

X [l + /(;i0)+/a20)], 

2.9^ 



d/, 



27r 



27r 



p{ho,h)pil2o, h)5 {po - ho - ho) 



X ed^i -^sl <P<h + h) 
X [l + /(;io) + /(/2o)], 
2.9^ 



p(27r) 



Zid^i 



hdh 



A-ho 
27r 



d^20 

27r 



X [(^10 - ^20)^] p{hoJi)p{ho, h)S {po - ho - ho) 
X &[\h-h\ <P<h+h] 
X [l + /(;io) + /a2o)] , 



Im n*p, = -7T 



29' 



hdh 



p{2n)- JO JO 
{ho — ho){li — I2) 



hdh 



d/10 
2tt 



dho 
2tt 



pihoJi)p{ho,h) 



X e - ^2! < P < ^1 + h) S {po - ho - ho) 
X [l + f{ho) + f{ho)] ■ 



(29a) 



(29b) 



(29c) 



(29d) 



With this set of equations for the self-energy and spectral-density components, we have established our self- 
consistently resummed approximation scheme. We solve this set of equations numerically. 



IV. FINITE WIDTH-EFFECTS ON THE p-MESON 



In order to gain insight into the physics employing the spectral modifications, we compare the result with the 
resummational approach of the previous section to perturbative one-loop calculations of the p- (cf. [13) ^^d n 
widths. 

As outlined in the last section, we work in the ^-functional approximation scheme in the i?j-gauge with ^ ^ 0, 
which is a choice guaranteeing decoupling of the Stiickelberg and Fadeev Popov ghost from the p vector meson field 
and fulfilling the transversality condition of the electromagnetic current (therefore respecting QED gauge invariance). 
Nonetheless as mentioned before this cannot cure the violation of the Ward Takahashi identities stemming from the 
violation of local gauge invariance by the ^-functional approach. That is why we also compare our method to the 
results obtained by the enforcement of the transversality condition by disregarding the temporal contributions of the 
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tensor (which is spoihng somewhat the ^-functional approach) which has been suggested and studied by |28l | in a 
similar model. We concentrate on the p-ir interaction. 

In order to gain information to what extend this influences the effects which we are interested in, we discuss the 
broadening effects in the full ^-functional approach (referred to as approach I) and in the transversality-enforcing 
resummation scheme (referred to as approach II). 

The phenomenological most important difference between the resummation approaches and the perturbative one is, 
that the perturbatively calculated self-energies of the p- meson show a threshold behavior: the self-energies for timelike 
momenta vanish if \/P^ < 2m7r, where P is the four- momentum of the meson relative to the rest frame of the medium. 
This leads to a characteristic threshold in the dilepton-spectra: only Icpton-pairs with invariant masses above 2m7r are 
emitted. Considering the self-energy components of the p-mcson in the self-consistent resummation scheme employed 
here shows that the finite width of the pions smears out this perturbative threshold behavior: the whole spectrum of 
invariant masses of the dileptons becomes accessible. Therefore this threshold ought to be considered as an artefact 
of the perturbative calculation. 

A. Self-energies and spectral densities 

We start our discussion with the four-dimensionally transverse self-energy components Imlla and Imllb of the p- 
meson. The ImHa-componcnt is spatial-transverse and the Imllij-component is spatial-longitudinal. These components 
are plotted as a function of energy for a fixed momentum p = 125MeV/c relative to the medium for three different 
temperature in fig. I^fleft and middle). 

In both self-consistent approaches the spatial transverse and spatial longitudinal imaginary parts of the self-energies 
show that there is no two-pion threshold. The perturbatively calculated imaginary parts of the self-energies can be 
compared to the self-consistent ones and vanish for < \/P^ = \/Po ~ < ^TOtt ■ 

The self-energy projectors are (up to possible poles and cuts) analytic functions in the complex energy-plane. 
Possible cuts on the real energy- axes corresponding to the imaginary part of the self-energy projections separate the 
upper and lower half-planes |27| . Because of the anti-symmetry of the imaginary-parts of the self-energy projectors 
these cuts extend in the self-consistent calculations (approach I and approach II) over the total axes. The retarded 
functions in the upper half plane are therefore totally separated from the advanced in the lower plane. This leads in 
both self-consistent approaches to a (static) dilepton-rate that has no two-pion-threshold behavior as discussed below. 

In approach I the Imllb-component shows another remarkable property. It has not only a change of sign at the 
light-cone (that means at the transition between time- and spacelike momenta for po ^ p): but also a singularity. The 
transversality-condition PTE^;^ means on that level that for = the component Inillb has to vanish. Indeed, that 
is not fulfilled in the <I>-derivablc approach (approach I), but is enforced in approach II. In approach II the self-energy 
component Imllt, remains a continous function of po even in the vincinity of the light-cone. 

Ean. ljAllbp shows that the energy component can have a 1/P^-pole that leads to a change of sign in Imllb at 
Pq ~ p 'AS obsvered in approach I. This is an artefact of approach I that is closely correlated to the violation of the 
transversality condition of the self-energy tensor II^j^, see also the discussion in section Hill Although transversality 
of the self-energy components is violated, the spectral-density tensor still remains four-transverse in the P^- 
gauge for ^ — > guaranteeing the conservation of the hadronic-electromagnetic current. Therefore approach I is 
consistent within the vector meson dominance approach with quantum electrodynamics, as already discussed above. 
The singularity of strength 1/P^ of the Imllij-component at the light-cone still leads to a continuous behavior of pf, 
at po = P- 

In figure |2| we plot the components of the spectral- function of the p-meson at a temperature T = 100 MeV as a 
function of energy po and momentum p as obtained within the self-consistent approach I. 

We plot the corresponding sp atial-tran sverse and spatial-longitudinal components of the spectral function of the 
p-meson depending on P^ = ^/p^ — p^ for a given momentum of p = 125 MeV/c for three different temperatures, 
the two self-consistent approaches and the perturbative calculation in figure^ 

The spatial longitudinal and transverse components differ i n bo th self-consistent calculations from each other only 
in the region about < -y/P^ < 500 MeV and are for about Vp^ > 500 MeV identical. With rising temperature this 
difference in the lower mass region is enhanced and the spectral components get in both self-consistent calculations 
considerably broadened. The perturbatively calculated spatial transverse and longitudinal spectral components (which 
are within the resultion of the figure almost identical) are considerably narrower and show no spectral strength below 
the two-pion-threshold. In figure we show the components of the spectral- function of the pion at a temperature of 
T = lOOMeV as a function of energy po and momentum p for the self-consistent calculation within approach I. Fig. El 
shows the tt spectral function for different T and a constant momentum of p = 125 MeV/c. It shows considerable 
broadening of the pion in both self-consistent approaches. Comparing to the perturbative calculation one oberserves 
apart from broadening the vanishing of the spectral function in the range between roughly 630 and 910 MeV. This 
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FIG. 2: The imaginary parts of the spatial-transverse and spatial- longitudinal self-energies Ha and 11;,, respectively, as a 
function of energy at a constant momentum of p = 125 MeV/c on the left hand side and in the middle. We show the results 
for a perturbative calculation and for the self-consistent calculations using the full ^-functional approach (approach I) and the 
transversality-enforcing approach (approach II). The imaginary part of the negative of the pion's self-energy is depicted on the 
right hand side. The perturbative self-energy of the pion is too small to be seen on the scale of the graph on the right hand 
side. 



is due to the fact that for two particles with mass mi > m2 the perturbative calculation leads to a vanishing width 
for mi — 1712 < ^J~P^ < mi + m2 - which corresponds to the two-pion threshold in the case of the p-meson and to the 
threshold behavior as mentioned here for the pions. This threshold behavior is no longer present in both self-consistent 
approaches. 
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FIG. 3: (Color online) The spatial-transverse and spatial- longitudinal spectral functions pa and pi,, respectively, of the rho 
meson as a function of energy and momentum on at a temperature of 100 MeV in the full ^-functional approach (approach I). 



B. Static thermal dilepton rate 

We proceed by calculating the static thermal dilepton emission rate in the model. For this purpose, we use ean. (|B5|l . 
where this quantity for the hadronic channel is given as a function of the invariant dilepton mass and as a functional 
of the p- meson spectral function based on VMD (see appendix IbI for details). We calculate the rate for different T 
(see figured for decaying p- mesons of momentum p = 125MeV/c relative to the medium. 

The dilepton emission rate increases with temperature. This enhancement is due to the Bose Einstein factor in 
eq. (|B5|) . The self-consistently calculated rates in approach I and II give rise to a larger yield in the low-mass part AI = 
\fp^ < 800 MeV of the spectrum than the perturbatively calculated rate which moreover has the two-pion threshold 
below 2m7r. In the present model, this leads to a considerable enhancement of the dilepton-rate especially in the \/P^ < 
700 MeV sector in comparison to the perturbative predictions. The self-consistently calculated rates of approach I 
exhibit even larger broadening effects for the p-meson than those of approach II. Despite this difference both self- 
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FIG. 4: The spatial-transverse and spatial- longitudinal spectral functions pa and pb, respectively, of the rho meson are given 
as a function of invariant mass at a constant momentum of p = 12f)MeV / c on the left hand side. We give the results for 
the perturbative calculation and for the self-consistent appoaches I and II (see text). The perturbative spatial-transverse and 
spatial-longitudinal components are within the resolution of the figure identical, that is why we used only one curve for both 
of them. 



-p in 1/MeV^ 




FIG. 5: (Color online)The spectral function of the pion as as a function of energy and momentum at a temperature of 100 
MeV as obtained within the self-consistent calculation applying approach I (see text). 



consistent approaches show clearly that in this part of the low mass dilepton region the coUisional broadening effects 
enhance the expected static rates considerably and outweigh the two-pion threshold of the perturbative approach. 
Above 850 MeV there are no important difference in comparison with the perturbative results. 
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FIG. 6: The spectral functions p of the pion is given as a function of invariant mass at a constant momentum of p — 125MeV/c 
on the left hand side. We give the resuhs for the perturbative calculation and for the self-consistent one. The width of the 
pion right hand side for the self-consistent and the perturbative calulation as a function of invariant mass. 




FIG. 7: The figure shows the static rate of dilepton (dielectron) production as a function of invariant mass. The decaying rho 
meson has a constant momentum of p — 125MeV/c. We give the self-consistent results as obtained within approach I and 
approach II (see text for explanation) and the perturbative ones. 



V. DILEPTON EMISSION AND FIREBALL EVOLUTION 



In order to demonstrate that the results of the previous section, especially the considerable coUisional broadening of 
the /5-meson leads to significant quantitative changes in the prediction of dilepton rates we model the fireball evolution 
and fold it with the rate calculations. 

We emphasize that we study a simplified model of the medium considering only the self-consistent finite width 
effects of the p-ir interaction at finite T where other effects (like mass shifts, chiral symmetry, and scattering off by 
baryons) have not been included. We do not suggest that these effects are unimportant but wish to demonstrate the 
self-consistent description of collisional broadening already implies siginificant modifications. 



A. The Fireball model 



In order to obtain the complete emitted spectrum we use a mode for the fireball which has been shown to reproduce 
a number of other important observables |l4| . Here, we briefly outline the framework. 

The fundamental model assumption is that the fireball matter is in local thermal (but not necessarily in chemical) 
equilibrium from an initial proper-time scale tq till a later breakup at time Tf. We chose volumes corresponding to a 
given proper time r for the calculation of thermodynamics. Assuming spatial homogeneity, temperature T, entropy 
density s, pressure p and chemical potentials fii as well as energy-density e become functions of r in such a system. 

The volume is assumed to be cylindrically symmetric around the z-axis of the beam. It is characterized by the 
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longitudinal extension L(t) and the transverse radius R{t): 

V{t) ^'!iL{t)R^{t). (30) 

In order to account for collective flow effects, we boost individual volume elements according to a position-dependent 
velocity field. For the transverse flow, we make the ansatz 

ilT{r,T)=r/R„ns{T)TlT'{r) (31) 

where Rrms{T) denotes the root mean square radius of the fireball at r and vi^^'^{t) the transverse rapidity at Rrms- 
For the longitudinal dynamics, we start with the experimentally measured width of the rapidity interval of observed 
hadrons 2r]j^°^* at breakup. From this, we compute the longitudinal velocity of the fireball front at kinetic freeze- 
out vj^°^^* . We do not require the initial expansion velocity Vq^°^* to coincide with vj^""^ but instead allow for a 
longitudinally accelerated expansion. This implies that during the evolution rj = r]s is not valid (with r]s the spacctime 
rapidity t]s = l/21n((t + z)/{t - z)). 

The requirement that the acceleration should be a function of r and that the system stays spatially homogeneous 
for all r determines the velocity field uniquely if the motion of the front is specified. We solve the resulting system 
of equations numerically 0|. We find that for not too large rapidities 77 < 4 and accelerations volume elements 
approximately fall on curves const. = \Jt'^ — z"^ and that the flow pattern can be approximated by a linear relationship 
between rapidity ry and spacctime rapidity 77^ as 77(778 ) = C'7s where C = 77'^''°"*/^/'^°"* ^-nd rj^"^""* is the rapidity of the 
cylinder front. In this case, the longitudinal extension can be found calculating the invariant volume V = J daf^u^ as 

L{r) « 2/^"^-«^^-^>;p*(-» (32) 

with r](^°"'*{T) the spacctime rapidity of the cylinder front. This is an approximate generalization of the boost-invariant 
relation L{t) = 2rjf'^°"*T which can be derived for non-accelerated motion. 



1. Parameters of the expansion 

In order to proceed, we have to specify the longitudinal acceleration az^r) (which in turn is used to calculate 
?7^''°"*(r) numerically), the initial front velocity Wq'^""* and the expansion pattern of the radius R{t) in proper time. 
For the acceleration we make the ansatz 

az = Cz ■ (33) 
e(r) 

which allows a soft point in the EoS where the ratio p/e gets small to influence the acceleration pattern. Cz and 
yfront model parameters governing the longitudinal expansion and fit to data. Hence the EoS translates into a 
temperature (and proper-time r) and determines the acceleration profile a = const. p{T) / e{T) . 

Since typically longitudinal expansion is characterized by larger velocities than transverse expansion, i.e. v[^°^* ^ 
^front ^ we treat the radial expansion non-relativistically. We assume that the radius of the cylinder can be written as 



i?(T) = Ro + CT I dr' I dr"^^ (34) 



TO 



The initial radius Rq is taken from overlap calculations. This leaves a parameter ct determining the strength of 
transverse acceleration which is also fit to data. The final parameter characterizing the expansion is its endpoint given 
by Tf, the breakup proper time of the system. 



2. Evolution of temperature 



We assume that entropy is conserved throughout the thermalized expansion phase. Therefore, we start by fixing the 
entropy per baryon from the number of produced particles per unit rapidity (see e.g. |34j). Calculating the number 
of participant baryons (see ^^^^ the total entropy Sq. The entropy density at a given proper time is then 

determined by s = Sq/V{t). 
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We describe the EoS in the partonic phase by a quasiparticle interpretation of lattiee data which has been shown 
to reproduce lattice results both at vanishing baryochemical potential /is and finite fxs [ssj . 

For the phase transition temperature, we choose Tc = 170 MeV based on lattice QCD computations at finite 
temperature for case of two light and one heavy quark flavour [s^ . We also note that no large latent heat is observed 
in lattice calculations for the transition and model the actual thermodynamics as a crossover rather than a sharp phase 
transition. Nevertheless, in the calculation we assume quarks and giuons as degrees of freedom above Tc and hadrons 
below to simplify computations. Since the time the system spends in the vicinity of the transition temperature is 
small compared with the total time for dilepton emission, any error we make by this assumption is bound to be small 
as soon as we consider the measured rates which represent an integral over the time evolution of the system folded 
with the emission rate. 

Since a computation of thermodynamic properties of a strongly interacting hadron gas close to Tc is a very difficult 
task, we follow a simplified approach instead: We calculate thermodynamic properties of a hadron gas at kinetic 
decoupling where interactions cease to be important. Here, we have reason to expect that an ideal gas will be a good 
description and obtain the EoS using an ideal resonance gas model. Using the framework of statistical hadronization 
|37j , we determine the overpopulation of pion phase space by pions from decays of heavy resonances created at Tc and 
include this contribution (which gives rise to a pion-chemical potential of order /^^r ~ 120 MeV into the calculation). 
We then choose a smooth interpolation between decopling temperature Tf and transition temperature Tc to the EoS 
obtained in the quasiparticle description. This is described in greater detail in [T^ . 

With the help of the EoS and s(t), we are now in a position to compute the parameters p(t), e(T),T{T) as well. 
Since the ratio p{T)/e{T) appear in the expansion parametrization, we have to solve the model self-consistently. 



3. Solving the fireball- evolution-model 

In order to adjust the model parameters, we compare with data on transverse momentum spectra and HBT 
correlation measurements. This is discussed in greater detail in [T3.l3^. 

By requiring R{Tf) ~ Rf and u^*^""* = v±f we can determine the model parameters ct and Tf. Cz is fixed 
by the requirement 77-^''°"* (tj) = 77^'^°"*. The remaining parameter Vq^°"^ now determines the volume (and hence 
temperature) at freeze-out and can be adjusted such that T{Tf) = Tf. 

The model for 5% central 158 AGeV Pb-Pb collisions at SPS is characterized by the following scales: Initial long, 
expansion velocity Wq'^""* = 0.54c, thermalization time tq = 1 fm/c, initial temperature Tq = 300 MeV, duration of 
the QGP phase tqcp — 6.5 fm/c, duration of the hadronic phase Thad — 8.5 fm/c, total lifetime Tf — tq — 15 fm/c, 
r.m.s radius at freeze-out Ry^'^ ~ 8.55 fm, transverse expansion velocity v±f ~ 0.57c. 

For the discussion of dileptons, we require the fireball evolution for other than 5% central collisions. In this case, we 
make use of simple scaling arguments based on the initial overlap geometry and the number of collision participants. 
For a detailed description, see [T^lsgj. 

In Il2|. it has been shown that this scenario is able to describe the measured spectrum of low mass dileptons, 
and in |37j it has been demonstrated that under the assumption of statistical hadronization at the phase transition 
temperature Tc, the measured multiplicities of hadron species can be reproduced. In [s^, the model has been shown 
to describe charmonium suppression correctly. None of these quantities is, however, very sensitive to the detailed 
choice of the equilibration time tq. Therefore, we have only considered the 'canonical' choice tq = 1 fm/c so far. The 
calculation of photon emission within the present framework provides the opportunity to test this assumption and to 
limit the choice of tq. In this has been investigated in some detail. Within the present framework, the limits 0.5 
fm/c < To < 3 fm/c could be found. Variations within these limits, however, do not affect the spectrum of dileptons 
with invariant mass below 1 GeV significantly. 



B. Dileptons from the QGP Phase 

In order to take into account the contributions from the phase of a QGP, we employ a quasiparticle picture |35j . 
The model treats quarks and giuons as massive thermal quasiparticles. The QCD dynamics is incorporated in the 
thermal masses of these 'effective' quarks and giuons. The masses in this model approach the hard thermal loop 
(HTL) results at high temperature. Around Tc a power law fall-off is assumed for the thermal masses, based on 
the conjecture that the phase transition is either weakly first order or second order. The thermodynamic quantities 
of the QGP are calculated in terms of two functions B{T) and C(T) introduced in the model, see ref. [mm , 
which account for the thermal vacuum energy and a phenomenological description of deconfinement. As long as the 
thermodynamically active degrees of freedom are quarks and giuons in a quasiparticle gas, the virtual photon couples 
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electromagnctically only to the thermally excited gq-statcs populating the plasma phase and converts into a lepton 
pair, for details see section 4.1 of This mechanism can be compared to the coupling of the virtual photon to the 
vector-meson in VMD in the hadronic phase. 

As the preferred reference frame of hot matter breaks Lorentz invariance, new partonic excitations such as longitu- 
dinal gluonic plasmons or helicity-flipped plasminos could be present. Their spectral strength remains exponentially 
suppressed for high temperatures and momenta which dominate macroscopic thermodynamic quantities (pressure, en- 
tropy and energy density) and their contribution is neglected in a quasiparticle description. However, for soft-dilepton 
production these plasmino modes lead to sharp, distinct structures in the static dilepton-emission spectra, so called 
van Hove singularities |42l | . These peaks should be roughly located at twice the thermal quasiparticle quark mass 
« 2mq (T) and would be smeared out by the fireball evolution, since the thermal quark masses drop with temperature. 
Furthermore, since in the quasiparticle model |35l | mq{T) is of order T (smeared out) van Hove singularities could only 
appear at low invariant masses and would be overwhelmed by the hadronic part of the dilcpton production. Therefore 
the effects of neglecting gluonic and (helicity-flipped) plasmino modes are expected to be neglible. 



C. Dilepton Emission 



Comparison with the measured dilepton spectrum is enabled by folding the differential rate dN / [d'^xd^q) = dR/d'^q 
with the space-time history of the collision. We want to compare the model predictions with the CERES/NA45 data 

Hi S taken in Pb-Au collisions at 158 AGeV (corresponding to a cm. energy of ^/s ^17 AGeV) and 40 AGeV 
{^/s ~ 8 AGeV). The CERES experiment is a fixed-target experiment. In the lab frame, the CERES detector covers 
the limited rapidity interval 77 = 2.1 — 2.65, i.e. Aij = 0.55. We integrate the calculated rates over the transverse 
momentum pT and average over 77, given that d'^p = Mpx dM drj dpx dO. The formula for the space-time- and 
p-integrated dilepton rate in the framework outlined above hence becomes 



d-'N 2nM 7, /• , ,,,, , ^, f, dNiT{T),M,r^',pT) , , , 

- ' ' dr] V{r],T{T)) J dpr pt ^cc(M, 77 ,pt), (35) 



dMdT] A77 



dr 



To 



where r/ is the freeze-out proper time of the collision. V^(7/,T(t)) describes the proper time evolution of volume 
elements moving at different rapidities and the function Acc{M,r],pT) accounts for the experimental acceptance cuts 
specific to the detector. At the CERES experiment, each electron/positron track is required to have a transverse 
momentum px > 0.2 GeV, to fall into the rapidity interval 2.1 < rj < 2.65 in the lab frame and to have a pair opening 
angle > 35 mrad. Finally, for comparison with the CERES data, the resulting rate is divided by dNch/di], the 
rapidity density of charged particles. 



D. Hadronic cocktail in the vacuum 



After freeze-out there are still hadrons present which decay after a finite time in vacuum and thus contribute to 
the dilepton yield. Below approximately 400 MeV Dalitz decays are the dominant decay mechanism in this regime. 
We use the experimental analysis of the CERES collaboration for SPS conditions. 

Most notably, remaining lo and (f> mesons contribute a characteristic sharp peak structure due to their vacuum 
decays. We calculate their contribution directly from the model. 

The direct decay of a vector meson V into lepton pairs is given via: 



dNy I , , f°° f M , n , f T -Tf 



Tf J y \ /vy^'o 

The temperatur Tf = T(Tf) and the volume Vf ~ ^(t/) of the the system are to be considered at freeze-out time 
Tf. After this time all medium effects vanish and the vacuum spectral function of the (virtual) photon Ry{M, T = 0) 
determines the rate of dilepton production. The momentum distribution of the vector mesons is given by the Bosc- 
Einstein contribution determined by the temperature at freeze-out. Due to the decay of the mesons and the absence of 
thermal recombination effects the number of vector mesons is decreasing exponentially as exp (— (t — T/)/(7(g)To')). 
Here, is the vacuum life time of the vector meson. The quantity 7(g) accounts for time dilatation effects on 
particles with finite three-momentum: 

( ^ 1 9" 
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FIG. 8: Dilepton invariant mass spectra normalized to dN^h/drj in units of (lOOMeV)"^ for the SPS-CERES/NA45 Pb(158 
AGeV)+Au (left side) and SPS-CERES/NA45 Pb(40 AGeV)+Au (right side) experiment. Shown are the data, the total rate. 
The theoretical spectra are calculated assuming perfect detector resolution. 



Integrating over proper time one obtains: 

The averaged four-dimentional space-time volume available after freeze-out is as one would expect: VfT^ . The three- 
momentum integral yields the particle density of virtual photons as a function of the invariant mass: n{M) = N{M) /V. 
Multiplying the factor Vf one obtains the total number of virtual photons decaying in the freeze-out stage. Weighting 
this with Rv{M) this gives the dN/dM distribution for the process in the freeze-out stage vector-meson 7* e^e~ . 



VI. DILEPTON INVARIANT MASS SPECTRA 



Once the time evolution is given in terms of the temperature T{t), baryon density p(t) and the knowledge of the 
dilepton rates from the hadronic medium at finite temperature and at freeze-out as well as the rates from QGP, all 
necessary ingredients for a quantitative calculation of the dilepton rates using eqn. H35|l are in place. We consider 
the SPS CERES/NA45 data taken at beam energies of 40AGeV and 158AGeV. The results for the dilepton invariant 
mass spectra for both experiments are shown in fig. |H1 

In the calculation we have assumed perfect detector resolution in invariant mass in order to demonstrate how 
the uj and (j) vacuum contributions would show although these structures in the spectrum cannot be resolved in the 
experimental analysis. 

Shown in Fig. IHlare the results of our model calculations as obtained in the two different self-consistent approaches. 
The solid line indicates the results in the full ^-functional approach (approach I), the dashed line the results in the 
approach where the transversality condition of the self-energy tensor is enforced (approach II). Both approaches lead 
to significant contributions for low invariant masses in the range of 200 MeV < Af < 700 MeV. This is due to the 
collisional broadening of the p- meson by the Tr-p interaction in the thermal medium which is taken into account by 
the self-consistent calculations. This principal finding remains in both self-consistent approaches the same. 

Comparing approach I and approach II one finds - beyond the confirmation of this in-medium collisional broadening 
- both results differ quantiatively in the 200MeV < M < SOOMeV region as already expected from the lower static 
dilepton rates of approach II in this mass region. These differences result from the two different possible methods to 
guarantee electromagnetic current conservation in the resummation schemes as discussed above which can ultimately 
be traced back to the lack of a feasible resummation scheme which guarantees local gauge invariance. We emphasize 
that this ambiguity in the results does not alter the principle finding, namely the importance of self-consistent 
collisional broadening effects induced by the ir-p interaction for an appropriate understanding of the observed spectra. 

In order to demonstrate that these broadening effects are quantitatively of the order of magnitude that is important 
for an experimental understanding, we compare them the results to the 40 AGeV and 158 AGeV data. We emphasize 
that many caveats have to be kept in mind comparing these results to data: we consider only the effects due to 
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coUisional broadening of the pions and the p mesons via their interaction in a simphfied model and have not taken 
into account effects like scattering off by baryons, the influence of chiral symmetry, or changes of the in-medium 
mass. Therefore we do not expect our simplified model to describe all details of the data. However, what we have 
shown with this comparison is that the self-consistent description of coUisional broadening, which has been neglected 
in the analysis so far has a strong influence on the theoretical description of the spectra in the dilepton yield in the 
low-invariant-mass 200 MeV < M < 750 MeV range. 

VII. CONCLUSION AND OUTLOOK 

We have investigated the in-medium properties of p-mesons in a self-consistent calculation taking into account 
finite-width effects due to the broadening induced by the n-p interaction at finite temperature. For this purpose we 
treat the p-mesons as a massive vector-mesons in a Stiickelbcrg formalism and resum the p — 7r-interaction considering 
finite-widths effects only. 

Special emphasis was given to the fact that this resummation in the Cornwall Jackiw Tomboulis formalism violates 
the Ward Takahashi identies connected to the Stiickelbcrg field and therefore gauge invariance. To circumvent these 
difhculties we have worked in a gauge which guarantees the transversality of the p-meson spectral-function in the 
resummation scheme. We also employed another approach first introduced by van Hees and Knoll in where 
they discussed self-consistent finite width effects on the p meson by pions in a similar framework. In that approach 
(approach II) one enforces the Ward Takahashi identities by disregarding the temporal components of the self-energy. 

We derived in detail self-consistently coupled Schwinger Dyson eqns. in the finite width approximation for the 
p-meson and the pions for both approaches and reformulated these equations into coupled integral-equations giving 
explicit expressions for the imaginary parts of the self-energies and the spectral-functions. 

Discussing the numerical solution of these equation gave physical insight into the broadening effects of the p-meson 
and the pion in a hot environment. 

This study demonstrates that p-mesons embedded in a pure pion gas exhibit substantial coUisional broadening, in 
contrast to perturbative approaches. This shows that a large broadening of the pions does not necessarily require a 
high density of baryons or chiral symmetry restoration. 

In order to demonstrate that these effects are quantitatively important and contribute to the understanding of the 
experimentally investigated dilepton spectra we made a comparison of the results with data obtained by the CERES 
NA49 collaboration using a parametrized fireball evolution model. 

We emphasize that this comparison does not aim at a description of all details of the data (where many important 
effects that we disregard in our model (as scattering off by baryons, chiral symmetry and resulting mass shifts) would 
have to be taken into account). 

We demonstrated that the investigated self-consistent finite-width effects on the p-meson are an in-medium modi- 
fication whose order of magnitude is not suppressed as compared to other effects. Therefore, they are relevant to the 
understanding of the experimentally oberserved enhanced dilepon yield (as compared to naive extrapolations of p-A 
collisions) in the low invariant mass region 200 MeV < M < 600 MeV in A-A collision. 

The phenomenological most important remaining questions might be how to handle and disentangle the additional 
effects of chiral symmetry restoration, interactions with baryons, and corresponding mass shifts additionally to the 
finite width effects on the p spectral function. 

Regarding the interaction with baryons, several studies already deal with the inclusion of medium modifications of 
vector mesons which arise from the modification of pionic modes modified through the nuclear medium at finite density 
and temperature Especially the last paper includes a treatment of the nucleon and the A(1232)-isobar 

resonance in a <I>-derivable approximation supplemented by a Migdal's short range correlation for the particle-hole 
excitation. 

From a theoretical point of view one might wish to include the handling of renormalization effects |l9l l45l prop- 
erly going beyond the finite width approximation and improve in the direction of feasible gauge-invariant resummation 
schemes which could also be interesting for the study of fundamental theories like QCD beyond HTL approximation. 
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APPENDIX A: PROJECTION OPERATORS FOR VECTOR PARTICLES 

In this section of the appendix, we introduce a projection operator^ formahsm that is suitable for the description 
of the propagating modes of vector particles in the medium. The notation is similar to |47j . In the presence of a 
heat bath as a special frame of reference, lorentz-invariance is broken by the medium. Therefore the spectral densities 
become functions of the energy po and the magnitude of the three-momentum |p| of the particles with four- momentum 
P — (P07P) relative to the medium. 

We introduce the projector E projecting onto the subspace parallel to P'^: 

E-'^^-p^. (Al) 
Then we introduce a vector which is orthogonal to P^: 



^ -J^ ) = (5^"^ - E^"^) /. , (A2) 

where = (0,p). Then N'^ = —p^p^/P^. The definition of the other three projectors is: 

_ N^ N-" ^ ^ p-^ + pi^ TV" , A'^" EE g''" - B"" - E'"' . (A3) 

The tensor A^'' projects onto the spatial transverse subspace of P^: 

A"" ^ ^ , A'^ = - {S'^ - p'p^) . (A4) 
The tensor B'^'^ projects onto the spatial longitudinal subspace of P^: 

B- = -^ , B- = -^ , B-.-^p^-. (A5) 

This can be written as: 

PA* 

B^^- = g^- + 5^1 - __ _ ffP , (A6) 

where for /i = or 1/ = the terms 5''^ and piPj don't contribute, elsewise i = /i and j = v. The projectors C and E 
have the following forms, respectively: 

COO ^ ^pW_ co.^ PPo{pl+p') ^. c-=2P^fjP, (A7) 

E"" = ^ , ^''-^f , E^^^^fp^ . (A8) 

The tensorial decomposition of self-energies, spectral functions, propagators etc. (represented here by the generic 
variable X) can be written as: 

X^'^(P) = X^(P) A^'' + X^{P) B''" + X^iP) + X^P) W . (A9) 

The scalar components of the projection X^, X'^ and X^ can be obtained by projecting X^^ onto these 
components. Using the abbreviations: 

X\P) EE i {S^^ -ffP) X^HP) , xHP) ^ P^X^\P)p, , (AlO) 



* N.b.: although the Uterature refers to the introduced operators as "projection operators", one has always to bear in mind that the do 
not project onto orthogonal subspaces. This becomes obvious considering eqns. iA13l . 
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one obtains the following equations: 



„2 



X^iP) 



P 
p2 

1 



2N2 p2 

1 

--pi 

X%P) = X^'^(P)E^, 



X^"{P) + 2 ^ X°XP) p, + ^ X\P) 



00, OA . Po±^x°\P)p, + X'{P) 
Pop 



= p2 [plX°'>{P) + 2p,pX°\P)p,+p'x'{P)] . 

If the tensor is four-dimensionally transverse (e.g. P^X^^" = 0) these equations simplify significantly: 

= -X\P) , 



X^iP) = -^X\P) 
Po 



X'>'\P) = ^X'{P) , 
Po 

X°\P)p, = -^X\P) 



Po 



X'={P) = 
X^P) = 



The inverse of a tensor X^'' (if it exists) is: 



X-\P) = X-i' -(P) A,,, + X-'- '^(P) B,,, + X-'- %P) C^, + X-'^ ^(P) E^, , 



where 



1 

T 



x-i-°(p) = 



,5 



and S = X'°X'= - P^N^ {X'^f. One checks exphcitly that indeed X'^'^X^^ = g"^. 



(Alia) 

(Allb) 
(Alio) 
(Alld) 

(Alle) 

(Al2a) 
(A12b) 

(Al2c) 

(A12d) 

(Al2e) 
(A12f) 

(A13) 



(A14) 



APPENDIX B: THERMAL DILEPTON PRODUCTION RATE IN A VMD MODEL 

In this section of the appendix, wc quote some important formulas for the dilcpton production rate via the decay of a 
thermal /o-meson in a VMD-model We use the formalism as in McLerran and Toimela |43| . Relating the electromagnetic 
current-current correlator 



W^M^^I d'xe^''-^e{x,)Y,{H I [J^(.T),J,nO)] I H) 

'' H 

to the photon propagator via 



-PEh 



leads to the dilcpton rate as a functional of the current-current correlator 



dR 



Sq 127r4g2 



2to? 
1 + ^ 



(Bl) 



(B2) 



(B3) 
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These equation are generally valid independent of the modelling of the hadronic electromagnetic current. 

In vector meson dominance (VMD) models is the retarded current-current correlator proportional to the retarded 
p- meson propagator (up to order 

W^^ = ^m^S^, (B4) 

The equation for the dilepton production rate from the hot hadronic medium as a functional of the p-meson 
spectral-function is now: 

dR 1 a^m^/ 2ml\ f 4^ ^, , 

Here, [dR/d^q) is the production rate of dileptons with a total invariant mass and the rest mass m^, of one lepton. 
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